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1. Introduction. 

Conformal representations of compact Riemann surfaces play an essential role in approximation theory; 
in particular, they serve to describe the so called ratio asymptotic of orthogonal polynomials on the real line. 

Let /i be a finite positive Borel measure whose compact support supp(/x) has infinitely many points and 
is contained in the real line. By qn{w) — XnW^^ + ■ ■ ■ + >co, with >f„ > 0, we denote the nth orthonormal 
polynonomial with respect to fi; that is 

J qn{t)qm{t)dfi{t) = Sn,m, 

where Sn.m denotes the Kronecker delta. 

It is well known and easy to verify that the sequence {q„},n > 0, satisfies a three term recurrence relation 

tqn{t) = a„+i(7„+i(t) + 6„g„(t) + a„9n-i(^), n>l, 

bn G R, a„ > 0. The following result establishes a close connection between ratio asymptotic of orthonormal 
polynomials and the limit behavior of the recurrence coefficients (see [6]). 

Theorem 1.1. The following assertions are equivalent: 

• lim a„ = a > 0, lim 6„ = b. 

n n 

• lim ^"'^'^^ — ^{w) := h \ ( ) — 1, uniformly on each compact subset o/C\ supp(/i). 

n qn[w) 2a V \ 2a J 
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// any of these conditions holds, then supp(/i) = [b~ 2a, 6 + 2a] U e, where e is at most a denumerahle set of 
isolated points m K. \ [5 — 2a, h + 2a] . 

The function z — ip{w) is the conformal representation of C \ [6 — 2a, b + 2a] onto the complement of the 
unit disk such that (p{oo) = oo,ip'{oo) > (C denotes the extended complex plane). It is a solution of the 
algebraic equation 

o w — b 

z + 1 = 0, 

a 

which lives on a two sheeted compact Riemann surface. The numbers a and b are determined by the first 
two coefficients of the Laurent expansion of ip at oo; in fact, (p{'w) = w/a — b/a + 0{l/w), w — s- oo. In turn, 
the algebraic equation which characterizes ip is easy to determine if the values a, b are given. Any measure 
for which supp(/^) = [b— 2a, b+ 2a] U e, where e is as described in the theorem, verifies this ratio asymptotic 
if /i' > a.e. on [b — 2a, b + 2a]. This is known as the Rakhmanov-Denisov theorem (see [3] and [H]). 

Recently, (see [1], [2], and [5j) results analogous to those stated above and the Denisov-Rakhmanov 
theorem were obtained for multiple orthogonal polynomials of Nikishin systems of m measures. In this 
case, the ratio asymptotic is described in terms of a conformal representation of an m + 1 sheeted compact 
Riemann surface onto the extended complex plane. The reader interested in these results can check the 
references given. We will not dwell into details and introduce directly the relevant Riemann surface. 

Let Afc, fc = 1, . . . , TO, be a system of bounded intervals of the real line such that n Afc+i = %,k = 
1, . . . , TO — 1. Consider the (to + l)-sheeted Riemann surface 



7^= U 7^fe, 

formed by the consecutively "glued" sheets 

7^o:=C\Al, 7^fc :=C\(AfeUAfc+i), /c = l,...,m-l, 7^™ = C\A„, 

where the upper and lower banks of the slits of two neighboring sheets are identified. In fact, this describes 
TZ as, a. branched cover of C. Denote by tt the corresponding projection. Notice that 7?. is a compact Riemann 
surface of genus 0; therefore, it is conformally equivalent to C. Fix / e {1, . . . ,to}. We are interested in the 
conformal representation t/"^'-*, of TZ onto C such that 

(1.1) il)^^\w)^w + 0{l), w^oo^^\ il^^^\w)^Ci/w + 0{l/w'^), w^oo^°'^, 

where w*-'-* denotes the point in the sheet I which projects onto e C. When m = 1 and Ai = [&— 2a, 6 + 2a], 

In view of our motivation coming from the ratio asymptotic of multiple orthogonal polynomials of Nikishin 
systems of measures, we would like to find an explicit expression for ■0'''' or, at least, an algebraic equation 
characterizing it, in terms of the end points of the intervals A^. This problem seems to be very difficult and 
unsolvable in all its generality, because it requires the solution of high order algebraic equations. For this 
reason, in the main part of the paper, we will restrict our attention to the case to = 2, where we have been 
able to achieve our goals (see Theorem 3.1 below). 

There is another way to look at the Riemann surface TZ. If -0 : 7?. ^ C is any conformal representation 
and (7 : C — *■ 7?. is its inverse, then 

(1.2) G:^Txog 

is a holomorphic function from the extended Riemann sphere onto itself and hence it is a rational function. 
Its degree (understood as the number of solutions of the equation G{z) = w for a generic w) equals the 
number of sheets of the branched cover n : TZ C, that is to + 1. The end points of the intervals A^ 
determine, in fact, the values of G at its critical points, where G'{z) — 0. Therefore, we are interested 
in determining a rational function G of a given degree by its critical values and, possibly, some additional 
topological characteristics of the branched cover it defines. 
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We remark that a related problem of classifying rational functions with real critical values has been studied 
recently by various authors, see [1^ and references therein. It also turns out that our problem is related to 
a problem of determining the coefficients of a special Schwarz-Christoffel mapping of the upper half-plane 
onto a special Riemann surface, see 21 below. 

The paper is divided as follows. In Section 2, we describe some general properties of the conformal 
representations of TZ onto C. Sections 3 contains the statement of the main results of the paper which are 
obtained for the case when m — 2. Their proofs are carried out in Sections 4 and 5. The final Section 6 is 
devoted to the development of a numerical algorithm for solving the system (Syst) which in turn allows one 
to calculate the functions ip^^^^ numerically (for m ~ 2). 

Acknowledgements. We would like to thank the referee for carefully reading the manuscript and for 
some helpful suggestions. We also express our gratitude to Daniel Estevez Sanchez, who helped us with the 
numerical calculations, using the MatLab system. His work was subsidized by a Grant for the Support to 
Study by the Autonomic Community of Madrid (2008). 



2. Some general properties 

Let V": ^ be any two conformal representations of Tl, onto C, then ^ o ■tp~^ is an automorphism of C. 
Consequently, there exist constants a, 5, c, d, with ad — cb ^ 0, such that 

*(w) - ^ ' 



cip{w) + d 



Therefore, if we have found a conformal representation -ip, any other 5* can be expressed explicitly in terms 
of ip and the values of at three distinct points (alternatively, in terms of the values of at two points and 
its first derivative at one of the two points). For example (see (|1.1|) ). 

V;W(u;) = C2V^''(u')/(V^''^(w)-V("(oo('=')), fc, Z € {1, . . . , m}, fc ^ Z, 

where C2 ^ is an appropriate constant so that the Laurent expansion of the right hand side at 00^*^^ has 
leading coefficient equal to 1. 

Let i/'fe \ ^ = 0,1,..., TO, denote the branches of corresponding to the different sheets of TZ. The 
function V^*-'' is the solution of the algebraic equation 

m m 

(2.3) n - '/'I- ^(^)) = ^"^^ + E "^■(^)^' = 0- 

k=0 j=0 

The coefficients aj{w) are the so called symmetric functions which are known to be entire functions on the 
complex plane and can be expressed in terms of the branches of ')/''■'■* through the Vieta relations. 

By the definition of ao(w) = (-1)"+^ IlS^ V'i' ' (w) has no singularity on C; therefore, according 
to Liouville's theorem it is constant. On the other hand, each aj{w),j = 1, . . . ,to, has a simple pole at 
cxd; consequently, it is a polynomial of first degree. In (|2.3p . we can solve for w and from what was said 
before we find that w = g{z) = p{z)/q{z) = {'4>^^^)~^{z)^ where p is a polynomial of degree to + 1 and q 
is a polynomial of degree to. The algebraic equation which defines TZ is irreducible, therefore, (p, (?) = 1; 
otherwise, ip^^'^ would satisfy a lower order algebraic equation in z. The poles of g (which are cxd and the 
zeros of q) are the points ip'"''\'^''''^)i fc = 0, . . . ,to. The zeros of g (which are the zeros of p) are the points 
^(0 (0^'^)), /c = 0, . . . , TO (without loss of generality we can assume that ^ U^pAfc). Since is single 
valued, all these points are distinct. In particular, the zeros and poles of g are simple. 

It is easy to verify that 



(2.4) %p^^\w) =tp^^){w), weTZ. 
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Indeed, let (t){w) :— tp^'-^w). 4> and V'^'-' have the same divisor; consequently, there exists a constant C such 
that (j) — Ctp^^K Comparing the leading coefficients of the Laurent expansion of these functions at oo^'^, we 
conclude that C — 1. 

In terms of the branches of , the symmetry formula (|2.4p indicates that for each fc = 0, 1, . . . , to, 
(2.5) :l\(AfcUAfe+i) ^1 

(Ap = Am+i = 0); therefore, the coefficients (in particular, the leading one) of the Laurent expansions at 
oo of the branches of "0^'^ are real numbers, and 



(2.6) Vi'^(i±) = V^i''(%) = ^i'|i(i±), t e Ak+i 

(the second equality is due to continuity). 

3. The results for to = 2. 

Here we consider the problem of explicitly determining the functions G and -0'-''' for the case of two intervals 
Ai,A2. Using an affine transformation, if necessary, we can assume that Ai = [—/_*,—!] and A2 = [1,A], 
where A, /i > 1. 

Theorem 3.1. The functions ^^'^'^ and V''-^^ can be computed as follows: 

^ H{a) + ^ ^' ^ 2H{a) 1 - G-i ' 

where G{z) := H{z)/H{a), 



(3.7) H{z) = h + z 

(3.8) /i=l(a + a)( 



Az Bz 



2aa 



1- z 1+z ' 
(«-a)'\ 



1 — aa 



(3.9) A = i(l- /?)(!- a)(l-a)(l- 6), B = i (1 + /?)(! + a)(l + a)(l + 6) , 

(3 and b are the solutions of the equation 

9 / N (a — a)^ 

(3.10) x^ + {a + a)x + \ ^-3 = 0, 

1 — aa 

verifying (3 < —1, b > 1, and a and a are the unique solutions of the algebraic system 

J 2 (a + a){?> — aa — a — q;)(3 — aa + a + a) + (A — /i)(a — aY = 
I (A + fif{a - a)6 = 4 (3 + aaf{l - aa)(2 + a + a){2 - a - a) , 



(Syst) 



verifying — l<a<a<l. 

Since H^^{w) is the solution of the cubic equation 

z^ -{w + A- B - h)z'^ -{l + A + B)z + w-h^O, 

the functions ■f/;^^-' and ■0^^^ can be computed explicitly if we know a and a. 

As follows from the proof, the rational function G, given by this theorem, is given alternatively by (|1.2p . 
where g is one of the conformal homcomorphisms of C onto TZ. As part of the proof of Theorem l3.1l we will 
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also obtain that G and the real numbers (3, a, a, b constitute a unique solution of the following system of 
relations 

'g=£, degP = 3, degQ = 2, 
— 1 and 1 are poles of G, 



(3.11) 



lim.^oo %i > 0, 



/3<-l<Q;<a<l<6, 
G'{(3) = G'{a) = G'{a) = G'{b) = 0, 
-^i^G{P), X = G{b), 
^-G(a) G(a) 1. 

Theorem 3.2. Consider the subsets o/R^ 

21 = {{a, a) : -1 < a < a < 1}, = {(A, fJ.) : \ > 1, > ^■ 

Then system (Syst) defines a one-to-one correspondence between these sets. Moreover, this mapping trans- 
forms the set {{a, a) : —1 < a < a < 1} onto the interior of fl, the sides {a — —1}, {a = 1} of the triangle 
21 onto the rays {A = 1}, {/i = 1} C fl, respectively, and the vertex (ao,ao) = (^l7l) o/ 21 onto the point 

(i,i)er!. 

If Ai and A2 have the same Euclidean length, then fj, — X, and we obtain a very simple result. 
Theorem 3.3. If \ = ji, then 

{l-a^f G-i + 1 



= i±|! (G-^ + 1) ^ 



2a3 



:2) 



8a3 G-i - 1 



where G{z) H{z)/H(a), 



H{z) ^ z- 



[X-a^fz 



(l + a2)(i_22) ' 

and a is the unique solution on the interval (0,1) of the biquartic equation 

+ (16A^ - 8)a^ + ISa"^ - 27 = . 
In this case, H~^{w) is the solution of the cubic equation 



n 2 '^^ ^ 

z — wz ^ — z + w = . 

1 + a2 



4. Proof of Theorem 13. II 

4.1. A geometric definition of G. Suppose intervals Ai = [— /i, — 1] and A2 — [1,A] are given. Then 
there exists a conformal homeomorphism g : C — > TZ such that G := n o g : C — > C is a rational function 
G — P/Q with real coefficients, degP = 3 and degQ = 2. In fact, g can be constructed as follows. Consider 
a conformal map g from the closed upper halfplane H onto the simply connected set 

n+ := {C e 7^o : 3^(0 < 0} y {C e 7^l : 37r(C) > o} |J {c G 7^2 : 5^(C) < o}. 

Notice that 7?.+ is "one half of TZ and coincides with the boundary of 7?,+ . We can choose g so that 

for some P<zi<a<a<Z2<bwe have g{{—oo,P]) = (—00, — ^]^^\ 9{[P, = (—00, — g([zi, a]) — 
[-l,oo)(0), g{[a,a]) = [-1, 1] W,^ ^([a, 22]) = (-oo,l](2), g{[z2,b]) = [A,oo)(2) and5([&,oo)) = [A,c^)«. 

By Schwarz's reflection principle, we can extend g to a conformal map (which we denote also by g) from 
C into TZ, which is symmetric with respect to the real line. Then G := tt o g : C — > C is holomorphic; 
therefore, it is a rational function G = P/Q, where P and Q have real coefficients, since G maps M into M. 
As G has degree 3 and has three poles (zi, Z2,oo), then a fortiori degP = 3 and degQ = 2. 
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G has two finite poles and four critical points on the real line. Taking G{diz + ^2) (with di > 0, ^2 G M) 
instead of G(z), if necessary, without loss of generality we can assume that the poles of G are —1 and 1. 
These normalizations define G uniquely. 

The critical points of G satisfy 

/3<-l<a<a<l<6. 

Moreover, G increases on (— oo,/3), (a, a), (&,oo), and decreases on (/?, — 1), (— 1,0;), (a, 1), (1,6). Therefore, 
G{a) < G{a); we also have = G(/3) < G(a) = -1 and 1 = G{a) < G{b) = A. 
It follows that G satisfies all the relations (|3.1ip . 

This way of constructing G shows that it is in fact a kind of Schwarz-Christoffel mapping from H onto 
the Riemann surface TV^ , which may be viewed as an unbounded polygon, whose angles are all equal to 27r. 

4.2. An auxiliary rational function H. In order to compute G, let us consider the rational function 
(4.12) H{z) := cG{z) ^ degP = 3, degQ = 2, 

Q{z) 

where c is a positive constant such that lim^^oo H{z)/z = 1. 

Then, the requirements (j3.1ip on G arc equivalent to the following conditions on the rational function H: 

H = ^, degPi =3, degQ = 2, 
— 1 and 1 are poles of H, 

11 = H{(3)/H{a), X = H{b)/H{a), H(a) = -H(a), 
l3<-l<a<a<l<b, 
H'{(3) = H'{a) = H'{a) = H'{b) = 0, 
Jim^^oo H{z)/z = 1. 

Since —1 and 1 are simple poles of i7, we get from the two last equations in ()4.13|) the equalities: 

^ {z-fi){z~a){z~a){z^b) ^ ^ , A , B 



(4.13) 



(z2_l)2 {z-lf {z + lf ■ 

There for 6 

rri . , ^ r (C-/3)(C-«)(C-a)(C-fe) AZ ^ BZ 

H[zj = h+ —2 —2 d( = h + z + 



((2-1)2 ^ 1-z 1 + z ' 

where h = H{0), so that ((3Jl) holds. 

4.3. The deduction of system (Syst). We are going to prove equations (|3.8p - (|3.10p and the validity of 
(Syst). In order to do that, we will express /3, 6, h, A and B in terms of a and a and write down equations 
on a, a. The choice of a and a as unknown variables seems to be the best one from the numerical point of 
view, because they are bounded (—1 < a < a < 1). 
From the equation for H'{z) we deduce that 

(4.14) (z -f3){z- a)(z - a)(z - b) = {z^ - if + A{z + if + B{z - 1)^ 
Replacing z = 1 and z = — 1 in (I4.14p . we obtain 

(4.15) 4A= (l-/3)(l-a)(l-a)(l-6), 4B = (1 + /3)(1 + a)(l + a)(l + 6) . 
Identifying the coefficients in ()4.14|) . it follows that 

(4.16) - (3-a~a-b = 0, 

(4.17) (ia + Pa + Pb + aa + ab + ab= -2 + A + B , 

(4.18) - Paa - I3ab - [3ab - aab = 2A - 2B , 

(4.19) I3aab=l + A + B. 
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Using ((4T71) and (fiTTO]) . we get 

(4.20) a{(3 + b) + a{fi + fe) + /36 + aa = -3 + (3aah. 
Relation (|4T6l) gives /3 + 6 = — (a + a) , and we deduce that 

(4.21) /3b= + + =-3+^"-")'. 

1 — aa 1 — aa 

Therefore, /? and 6 can be obtained from a and a as the solutions of the equation 

4.22) x^ + ia + a)x + \ ^-3 = 0. 

1 — aa 

Since 

H{z) = h + z + - + -— , 

1 — z 1 + z 

the condition H{a) = —H{a) means that 

(4.23) Ah + Aa + a{l - /3)(1 - a)(l - 6) + a(l +/?)(! + a)(l + 6) 

= -4/i - 4a - a(l -/?)(!- a)(l - b) ~ a(l + /5)(1 + a)(l + b). 
Using that f3 + b = — (a + a), we express h in terms of a and a: 

(4.24) (a + a)(2aa~3- (3b) = (a + a) f 2aa - 

V 1 — aa / 

We are going to make use of the equations 

iXH{a)^H{b), 
\fiH{a) = -H{f3). 

Since H{P) < < H{b), these equations are equivalent to the following ones: 

{\H{P)+ ^iH{b)^0, 
\ {X + fi)Hia) ^H{b)~H{/3). 



(4.25) 



(4.26) 
We have 



2H{a) = 2/1 + 2a + I (1 - /3)(1 - «)(! - 5) + ^ (1 + /3)(1 + a)(l + b) 



^2h + a(3 - Q;(a + a) + bj3) , 
2F(6) = 2/i + b[3{a + a) + 6(3 + aa) , 
2H{(3) = 2h + b[3{a + a) + /3(3 + aa) . 
So we get from (|4.26p that 

f(A + /i)(2/i + &/3(a + ci')) + (A/3 + H(3 + aa) =0, 
^ ■ I (A + /i)(2/i + a(3 - a{a + a) + b(3)) (5 - /3)(3 + aa) . 

Since b and /? are the solutions of the equation + {a + a)x + b(3 — Q , we have 



-(a + a) + J (a + - 46/3 , ^ -(a + a) - J I a + a)2 - 4&/3 

= and p = 

2 ' 2 

From this, we deduce that 

b- (3^ v/(a + a)2 - 46/3 , 

and 

A + /i fi — X 



A/3 + /i6 = ^ {a + a) + ^(a + a)^ " 46/3 
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aa) — , 



(one could, in fact, express here &/3 in terms of a and a, see (|4.21|) ). These last two equations and (|4.2H) 

imply that (j4.27l) is equivalent to: 

(4.28) 

(A + n)Uh + (2 - 6)(a + a)) + ( - (A + fi){a + a) + {^l ^ A)^(a + a)2 - 46/3 )(3 - 

(A + /^)(2/i + a(-a(a + a) + if3^)) = (3 + aa)V(a + a)^ - 46/3 . 

After replacing Ah by its value, given in (|4.24p . we get 

/ (a — a)^ \, , , , , (a — a)^\ 

4/i+ 2^ 6 (a + a)- (a + a)(3 + aa) = (a + a)( - 9 + aa + —] 

\ 1 — aa / V 1 — aa/ 

and 

(4.29) 2g(a)^2fe + a(-a(a + a)+ = ~ ■ 

V 1 — aa / 2(1 — aa) 

Then, (|4.28p is equivalent to 



(4.30) 



(A + Ai)(a + a)(-9 + aa+if^) + (/i - A)v/(a + a)2 - 46/3(3 + aa) =0, 



(4.31 



+ A^) = 2 (3 + aa)V(a + a)2-46/3 . 

Substituting the second equation in (|4.30p into the first one, we obtain an equivalent system: 

2(a + a)(-9 + aa+(H?)+(A.~A)(f^=0, 
+ = 2 (3 + aa)V(a + a)2 -46/3 . 



The first equation in (|4.3H) is equivalent to the first equation in (Syst). Since — 1 < a < a < 1, both terms 
in the second equation in (|4.3ip arc positive. Therefore, this equation is equivalent to 

+ A^)'7F^ = 4 (3 + aa)' ({a + a)^ + 12 - 4 if^) , 
(1 — aa)^ V 1 — aa / 

and hence to the second equation in (Syst). 

We conclude that the functions G and H are necessarily given by the formulas in Theorem 13.11 where a 
and a are some solution of (Syst). The uniqueness of the solution will be checked in the next subsection. 

4.4. The reverse arguments. First we need some technical lemmas. 

Lemma 4.1. Let ti,t2 be the solutions of the real equation t' +pt + q = 0, with 3?ii < 3?t2- Then ti < ~1 
and t2 > 1 if and only if I + \p\ < —q. 

Proof. It is obvious because the condition ii < — 1 and ^2 > 1 is equivalent to /(—I) < and /(I) < 0, 
where f{t) :== t'^ + pt + q. □ 

Lemma 4.2. Let a, a be any real numbers such that —1 < a < a < 1. Then the solutions /3, 6 of the 
equation 

2 / N {o- — a)^ 

x' + (a + a)x + 3 = 

1 — aa 

are real. Assuming that (3 <b, one has /3 < — 1, 6 > 1. 
Proof. By Lemma HTTl it is sufficient to show that 

(a - a)2 

1 + |a + a| < 3 - — , 

1 — aa 

whenever — l<a<a<l. By symmetry, it suffices to verify that 

Ty(a,a) :=a + a+ <2, if |a| < a < 1. 

1 — aa 
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But this is a direct consequence of 



I — aa 



□ 



Lemma 4.3. Suppose X > 1, fi > 1 are arbitrary real numbers. Let G be the unique rational Junction that 
satisfies (j3.1ip . and define H from (|4.12p . 

(i) There is at least one solution {a, a) to the system (Syst) such that —1 < a < a < 1. 

(ii) Given any such solution, define [3 < —1, b > I from (|3.10|) . A, B, h from (|3.9p . (|3.8[) . and then 
define the corresponding function H from (|3.7p . Set G{z) — H{z)/H{a). Then G = G , H ^ H. 



Proof. Part (i) has been already checked in subsection [ 

Now assume that a, a, b, (3, H, G has been constructed as in (ii). We show that 

It suffices to check equations (|4.16l) - (|4.19p . 

We have (|4.16p from the definition of /? and b. This definition also gives (I4.2ip , which implies (I4.20p . This 
last equation is the difference of (|4.17p and (|4.19p . The definition of A and B gives 

2A + 2B = 1+ f3a + f3a + f3b + aa + ab + ab + /3aab , 

and this equation is the sum of (|4.17p and (|4.19p . Therefore, we can also obtain the equations (|4.17p and 
(|4?T9l) . The definition of A and B also gives (|4T8l) : 

-2A + 2B^l3 + a + a + b + l3aa + (3ab + (}ab + aab 

— (3aa + (3ab + (3ab + aab . 

Consequently, we have proved (|4.32p . 
The definition of h and formulas 

(a - 

P + b = ~a-a, 13b = ^ - 3 , 

1 — aa 

give (|4.23p . which is equivalent to H{a) = —H{a). 

We can reverse the arguments given in Section UT] to show implications (Syst) (I4.30p ^ (|4.26p . We 
also see that (|4.29p holds, which, together with the second equation in system (|4.26p . implies that H{a) ^ 0, 
H{b) ^ H{P). Therefore, we can deduce from ((T^ . Hence all the conditions hold. 

Now, we can assert that H increases on (— oo,/3), {a, a), (b,oo), and decreases on (/3, — 1), (— 1,q;), (a. 1), 
(1, b). Since H increases on {a, a), we have H{a) < H{a), but H{a) = —H{a), so we deduce that H{a) > 0. 

Then iJ'(/3) = H'{a) = H'{a) = H'{b) = 0, H{W} = R and lim,,^oo H{z)/z = 1. Furthermore, H has 
just three poles; namely —1, 1 and oo, and they are simple. 

The function G{z) = H{z)/H{a) verifies G(/3) = -n, G{a) -1, G(a) = 1 and G(6) = A. We also have 
G'(/3) = G'{a) = G'{a) = G'{b) = and G(l) = S. 

Since H{a) > 0, the function G increases on (— cx),/?), (a, a), (&, oo), and decreases on (/3, — 1), (— l,a), 
(a,l), (1,6). 

G is a rational function of degree 3; therefore, we conclude that G gives rise to a conformal map g of C 
onto a Riemann surface TZ* with three sheets (7?.*, 7?,*, T^-j), so that G = n o where tt : TZ* ^ C is the 
canonical projection. The set of branch points of TZ* is {G(/3), G(a), G(a), G{b)} = {—fJ., —1, 1, A}. All these 
branch points have order two, because G' has simple zeros on f3, a, a, b. 

Since TZ* = 7?.* U 7?.5; U TZ^ is connected and the branch points have order two, there is some TZ* (for 
instance, TZl) with two cuts and TZ*,TZ2 have just one cut. By the monotonicity properties of G, we can 
deduce that G(oo) = oo(i), G(-l) = oo^°\ G(l) = 00(2), and that G(H) = (7^*)+, where (7^*)+ is the union 
of one half of each 7^* (j = 0, 1, 2). We also have G(l) = d{TZ*)+. Every point of 1\ ([-/i, -1] U [1, A]) has 



10 



G. LOPEZ, D. PESTANA, J. M. RODRIGUEZ AND D. YAKUBOVICH 



3 preimages by G in M, and every point of (— /i, —1) U (1, A) has one preimage by G in M. Since (7?.*)+ is the 
union of one half of each TZ* and it is connected, the cuts are (— /i, —1) U (1, A). Consequently, TZ* — TZ. 

Therefore, G^^ o G is a conformal map from C onto C, and it must be a Mobius transformation. Since 
G(oo) = G{oo) = oo(i), G(-l) = G(-l) =j>o(o), G(l) = G(l) = oo^^), the Mobius map G'^ o G fixes the 
points 00, —1 and 1. It follows that G^^ o G is the identity map, and therefore G = G. □ 

Corollary 4.4. The solution (a, a) to system (Syst) with the properties — l<Q;<a<l is unique. 

Proof. Indeed, as shown above, any solution (a, a) of this system gives rise to a function G, which coincides 
with G. Denote the parameters that correspond to G by a, a, etc. Then a = a and a = a (because they are 
the only two critical points of G and G on the interval (—1, 1)). □ 

Since G{z) — H{z)/H{a), we obtain that G^'^{w) = H-^{H{a)w). We have 

rr, ^ + - + ^(1 + ■Z)^ + B{1 - Z)Z 

w = H{z) = ^— 5 , 

or equivalently, 

(4.33) z^-{w + A-B- h)z^ -{l + A + B)z + w-h = Q. 

Observe that (|4.33p allows one to obtain an explicit expression for H^^ (and then for ■0'-^'' and V'''^^); once 
we have solved (Syst). 

We compute now -0'^^ and -0'^^ in terms of G. Recall that -0'^^ and 0^^^ are determined by 

i;^^\w) =1(7 + 0(1), w^oo^^\ '0(i)(oo(°)) =0, 
ij^'^\w) =1(7 + 0(1), it;^oo(2), -0(2)1-^(0)^ 

Note that G-i(oo(")) = -1, G-^{oo^^^) = oo and G-^{oo^^^) = 1. Hence, V^^^ = ci(l + G^^) for some 
constant ci . In order to calculate the constant ci , notice that 

i7(z) = z + 0(l), z^oo, 
H~\w) = w + 0{l), w^oo^^\ 
G-\w) ^ H{a)w + 0{1), w^oo^^\ 
w + 0{l)=ciH{a)w + 0{l), w^oo^^\ 

and then 

_ 1 

We also have t/j^^'^ = C2(l + G~^)/(l — G~^) for some constant C2. In order to compute the constant C2, 
let us expand the left hand side of (|4.33|) in powers of (z — 1). One gets that (|4.33|) is equivalent to 

{z-lf + {3-w-A + B + h){z-l)'^ + {2-2w-3A + B + 2h){z - 1) - = , 

(z-l)2 {2,-w- A + B + h){z-l) 2-2w-3A + B + 2h 2A 



WW w u]{z — 1) 

Since ^^"^(00*^^^) = 1, we have that 



= 0. 



2A 2A 
lim — ■ — — -- = lim 



„^oo(2) w(l - i/"l(w)) ui^oo(2) w(l - z) 



or, equivalently, 



2 2 2H(a) 2H(a) 
hm — „ ., , — -- = lim — : 7- — — — -- = lim — 7- — -- = : — 

~+oc>(2) w{l - G-^{w)) «,^co(2) w{l - H-^{H{a)w)) w^ca(^) w{l - H-^{w)) A 
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Recall that 



i^^^Hw) l + G-\w) 2 2H{a) 

1 = iim — C2 iim — — -- = C2 iim — — -- = C2 



UJ — *oov- 



d(2) w u.^oo(2) w(1 - G"1(w)) u>^oo(2) w(1 - G"i(it;)) yl 

and this gives C2 = A/2H{a). This finishes the proof of Theorem 13. II 



5. Proofs of Theorems 13.21 and 13.31 
Proof of Theorem \3.2l We associate with any point (A, /i) G the point {u, v), where 
(5.34) u = X-fi, v = {X + nf. 

In terms of u, v, (Syst) may be rewritten as 

2 {a + a){3 — aa — a — a)(3 — aa + a + a) 



(5.35) 



(a — a)'^ 

(3 + aaf{l - aa){2 + a + a)(2 ~ a - a) 



(a-a)6 

One can check that these equations imply thalQ 

^^ (a2-l)(a2 + 2aa-3)3 (a^ _ i)(a2 + 2a« - 3)^ 

(5.36) V — [2 — u) =16 7^ , w — (2 + u) =16- 



Notice that -|- 2aa - 3 < for aU (a, a) e (-1, 1) x (-1, 1). It also follows that + 2aa - 3 < for all 
(a, a) e [-1,1] X [-1,1]. 

Take any point (a, a) e 21. By ()5.36p . the corresponding pair {u, v), given by ()5.35p . satisfies w > (2 — u)^, 
w > (2 + u)2. Hence w > 0, and y/v>2-u, y^>2 + u. Put 

(5.37) A=i(u + \Ay), /i=i(y^-M); 

then A, /z > 1. By examining (|5.34p . one gets that the pair (A,/i) is the unique solution of (Syst) that 
belongs to fl. By repeating the same arguments with strict inequalities, one also sees that A > 1, /i > 1 
whenever (a, a) £ int 2t. The assertions of the Theorem about the images in of the sets {a = —1}, {a = 1}, 
and {a = — 1, a= 1} follow easily from (I5.36p . 

It was proved already in the previous section that for any (A,/i) G intf2, (Syst) has a unique solution 
(a, a) e int 21. This finishes the proof of the Theorem. □ 

Proof of Theorem \3.3[ Assume that the intervals [— /i, —1] and [1, A] are symmetric, i.e. A = /i. Then, taking 
into account that (3 — aa — a — a)(3 — aa + a + a) > 0, we deduce from the first equation in (Syst) that 
a — —a. Hence (14.241) gives ft. = 0. By (|4.22p . it follows that /3 and b are the solutions of the equation 

4a2 3 - a2 



= 3- 



Consequently, /? = —b. Hence, by (|4.15|) . B — A = —777;—; — • Therefore, 



(l-a2)2 
2(1 + a2) 



(l-a^)^^ 2a3 

H{z) = z- H(a = — 

(1 + a^)(l — z"') 1 + 

and 



2a3 ' " 8a3 



^We found these formulas by using the package Maple of symbolic calculations. 
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Consequently, H is an odd function, and the second equation in (Syst) yields 

IGX^a^ = (1 + a2)(3 - = 27- 18a^ + Sa^ - a^. 
We have obtained that a is solution of 
(5.38) a^ + (16A2-8)a*^ + 18a'* -27 = 0. 

In this simpler case, we can check directly that there is exactly one solution of (|5.38p on the interval (0, 1). 
In fact, the function u{t) = t"^ + (IGA^ - 8)t^ + ISi^ - 27 verifies u(0) = -27 < 0,u(l) = IGA^ - 16 > 0, 
and u'it) = At^ + 3(16A2 - 8)t'^ + S6t > 36i > for every t e (0, 1). Furthermore, since ((OHl) is a biquartic 
equation, it can be solved by radicals. □ 

6. Numerical resolution of the nonlinear system. 

The purpose of this section is to show that the problem of solving numerically the algebraic system (Syst) 
is a simple task, both from the theoretical and practical points of view. 

Let us use the vector notation: A = {a, a) G int 21 and L = (A, /i) G int fl. By Theorem l3.21 (Syst) defines 
uniquely a map 

F : L. 

We suggest to use a continuation method for solving this equation numerically, which is a variation of 
methods described, for instance, in [4], Section 4.4.2. Notice that these methods go back to the treatise j8] 
by H. Poincare on celestial mechanics. 

Recall that the Newton method for solving equation F{A) = L with the starting point Ast consists in the 
following (see, for instance, |4j. Section 4.2). Fix a small parameter a > and define iteratively A^^^ = Ast', 
^(n+i) ^ ^(n) _ n > 0. This iterative process finishes when |F(A(")) - L| < ct, and the 

last point Afin := A*^") is taken for an approximate solution of equation F{A) = L. 

In the continuation scheme, we apply the Newton algorithm several times. Given a point i» G intfi, we 
wish to find a good approximation for a solution A^ G int 21 of the equation F{A^.) = L,. The method goes 
as follows. 

Part 1. Choose an initial approximation Aq of the solution. We do it solving the equation F{Aq) = Lq, 
where Lq = (A* + A* + /i*)/2 is a symmetric vector, by applying Theorem 13.31 

Part 2. Choose a large integer n > 0, and divide the interval [Lq,L^,] C into n equal subintervals by 
division points Lq, Li, . . . , Ln = i*. So, we set L^ = [{n — k)Lo + kL^^ /n. 

Part 3. Let the parameter cr > be fixed. The calculation is performed in n steps, applying the Newton 
method n times. 

On the fcth step (1 < fc < n), we find an approximate solution A^'^^^ of the equation F{A) = Lk by 
running the above-described Newton method with the starting point A^^^ = ^/^n^^' which is available from 
the (A; — l)th step, li k = 1, then we put A^^t = ^O: where Aq is the value found from Part 1. 

Once all n steps of Part 3 are performed successfully, one takes ^y"^ for an approximation of the solution 

A* of the equation F{A^) = L*. One has, in fact: -F(Aj"^) « L„ = L*. 

The applicability of this method for our concrete function F is justified by the following theorem. We 
will formulate it in a more general setting. 

Theorem 6.1. Let 2t and Vl be open subsets o/R'*, and suppose that 21 is connected and H, is convex. Let 
F be a homeomorphism from 2t to such that det F' ^ Q in 21. Then the above continuation scheme for 
the Newton method is numerically feasible in the following sense. 

Let a point G intf2 and a starting point Lq be given. Let G 21 6e the (unique) solution of the 
equation F{A) = L,. Then for any e > there exist a ctq > and an integer Nq > 1 such that for any 
n > Nq and any a < aQ, the Newton method stops on each of n steps of Part 3 of the algorithm, and the 
approximate solution ^i"^ — ^/"i(o') obtained satisfies — A^\ < e. 
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Because of the topological assumption on F, the situation is very simple. The feasibility of the method 
is closely related with the continuation property for F, see Section 5.3. 

It seems that this theorem is a variation of classical results, known to specialists. Since the authors were 
unable to find an exact reference, a sketch of the proof is included. 

Proof of Theorem \6.1\ Denote by B{q,r) the open ball centered at a point g G M'' of radius r. We apply 
Theorem 4.10 from [4]. It gives sufficient quantitative conditions for convergence of the Newton method, 
which imply, in particular, the following. Suppose A € 21, and let L = F{A) G CI. Then, by the assumption, 
detF'{A) 0. It follows that there exists S > with the following property. 

(*) Given any e > 0, there exists cto > such that to any a < there corresponds an integer 
M(cr) satisfying the following: the Newton method for solving equation F{A) — L with any starting 
point Ast S B{A,S) stops after at most Af(cr) steps, and the approximate solution Afin{a) obtained fulfils 
\Af,r,{c7) - A\ < e. 

Next, the pre-image F~^{[LQ,L^,]) is compact, hence the infimum of | detF'j on this set is positive. It 
follows that the above condition (*) takes place uniformly for all pairs of points {A, L) G 21 x such that 
F{A) = £ and L e [Lq.L^]. 

Fix Lq, L*, and fix some 5 such that (*) holds for all pairs {A^L) as above. Take any e > 0, and let 
us prove that the conclusions of the theorem hold for this e. We can assume that e < 5/2. Let us find 
the corresponding ctq such that (*) holds. There exists a positive p such that \A' — A"\ < e whenever 
F{A'),F{A") e [Lq,L4 and \F{A') - F{A")\ < p. Let iVo be any integer such that hNo > - Lo\. We 
claim that the assertion of the Theorem holds for these ctq and ^o- 

Indeed, take any n > Nq, and put (as above) Lk = [{n — k)Lo + kL^'j/n, and Ak — F^°-{Lk). Then 
\Ak — Ak+i \ < e for all k. Consider the properties: 



(r), \A^^A,\<S; (A), \Afl~A,\<e. 



Since A^j^ — Aq, (T)^ holds. Next, for any k, (T)k implies (A)^, due to (*), and (A)^ imphes (r)fe+i, 

- Afe+il = |4^i - Ak+i\ < \A%^ ~A,\ + \A, - Afe+il <2e<d. 
By induction, we obtain that (r)fc and (A)/c hold for all k. In particular. 



\A^;l-A.,\^\4t{a)-A.\<s. 



□ 



Let us return to our particular function F, defined by system (Syst). It follows from Theorem 13.21 and 
explicit formulas (I5.35p . (I5.37P that F : int2t ^ intfi is a C°° smooth homeomorphism. Put 2t = intSl 
and n = intri. In order to prove that the above Theorem 16.11 applies to F, it only remains to check that 
detF' does not vanish in int2l. This follows from an explicit calculation, using formulas (|5.35p . (|5.37p . that 
represent as a composition map. One has 

detF' = det§^=det|^.detf^. 

o[a,a) o[u,v) 0(a,a) 

Then, det = l/(4wi/2) > 0, because V = v{a,a) > in int2l. At last, (|5.35p implies that 

det = -2^ (a - a)-i0(a2 + 2aa - 3)^(3 + aa)^{a^ + 2aa - 3)^ 

d[a, a) 

(we have applied the Maple package here). As we noted in the proof of Theorem 13. 2[ + 2aa — 3 < and 
+ 2aa - 3 < for all (a, a) G int 2t. Therefore det |^ < and det F'{A) < for aU A = (a, a) £ int 21. 
In the next table, we reproduce some numerical results of the implementation of the above method. 



14 



G. LOPEZ, D. PESTANA, J. M. RODRfCUEZ AND D. YAKUBOVICH 

Table 1. Numerical results 



A 




/3 


a 


a 


b 


1.01 


1.10 


-1.02433457 


-0.97566543 


0.99756619 


1.00243381 


1.01 


1.50 


-1.11120778 


-0.88879220 


0.99776894 


1.00223104 


1.01 


2.00 


-1.20162987 


-0.79837009 


0.99796298 


1.00203698 


1.01 


5.00 


-1.53512664 


-0.46487320 


0.99855894 


1.00144090 


1.01 


10.00 


-1.80340393 


-0.19659585 


0.99893546 


1.00106432 


1.01 


20.00 


-2.05587223 


0.05587246 


0.99922942 


1.00077035 


1.01 


50.00 


-2.33815403 


0.33815423 


0.99950547 


1.00049433 


1.01 


100.00 


-2.50602945 


0.50602961 


0.99964857 


1.00035127 


1.10 


1.50 


-1.10904414 


-0.89095444 


0.97812791 


1.02187067 


1.10 


2.00 


-1.19803860 


-0.80195744 


0.97999986 


1.01999618 


1.10 


5.00 


-1.52833535 


-0.47164952 


0.98579364 


1.01419123 


1.10 


10.00 


-1.79584373 


-0.20413551 


0.98948341 


1.01049582 


1.10 


20.00 


-2.04866217 


0.04868414 


0.99237698 


1.00760105 


1.10 


50.00 


-2.33227827 


0.33229733 


0.99510249 


1.00487845 


1.10 


100.00 


-2.50132594 


0.50134152 


0.99651797 


1.00346645 


1.50 


2.00 


-1.18411168 


-0.81581475 


0.90730105 


1.09262537 


1.50 


5.00 


-1.50111605 


-0.49858460 


0.93306652 


1.06663413 


1.50 


10.00 


-1.76499059 


-0.23458419 


0.95001591 


1.04955887 


1.50 


20.00 


-2.01887210 


0.01933234 


0.96355662 


1.03598315 


1.50 


50.00 


-2.30775501 


0.30816158 


0.97647774 


1.02311569 


1.50 


100.00 


-2.48160429 


0.48193938 


0.98324158 


1.01642333 


2.00 


5.00 


-1.47250352 


-0.52657019 


0.87477665 


1.12429705 


2.00 


10.00 


-1.73155665 


-0.26707834 


0.90556526 


1.09306973 


2.00 


20.00 


-1.98589261 


-0.01259181 


0.93068264 


1.06780179 


2.00 


50.00 


-2.28012155 


0.28148839 


0.95501596 


1.04361720 


2.00 


100.00 


-2.45919974 


0.46033681 


0.96787312 


1.03098981 


5.00 


10.00 


-1.59758021 


-0.39392722 


0.70769711 


1.28381032 


5.00 


20.00 


-1.84539761 


-0.14406679 


0.77852683 


1.21093757 


5.00 


50.00 


-2.15584295 


0.16630302 


0.85234419 


1.13719574 


5.00 


100.00 


-2.35580726 


0.36490617 


0.89323258 


1.09766851 


10.00 


20.00 


-1.70949675 


-0.26791189 


0.61012964 


1.36727900 


10.00 


50.00 


-2.02314506 


0.04811583 


0.73131282 


1.24371641 


10.00 


100.00 


-2.239<S3()<S<S 


0.20275325 


0.8025()032 


1.17452331 


20.00 


50.00 


-1.86312347 


-0.09255784 


0.56794501 


1.38773629 


20.00 


100.00 


-2.09037160 


0.13418100 


0.67474116 


1.28144945 


50.00 


100.00 


-1.85498816 


-0.06912394 


0.44519539 


1.47891671 
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